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Dense neutrino gases exhibit collective oscillations where "self-maintained coherence" is a char- 
acteristic feature, i.e., neutrinos of different energies oscillate with the same frequency. In a non- 
isotropic gas, however, the flux term of the neutrino-neutrino interaction has the opposite effect of 
causing kinematical decoherence of neutrinos propagating in different directions, an effect that is 
at the origin of the "multi-angle behavior" of neutrinos streaming off a supernova core. We cast 
the equations of motion in a form where the role of the flux term is manifest. We study in detail 
the symmetric case of equal neutrino and antineutrino densities where the evolution consists of 
collective pair conversions ("bipolar oscillations"). A gas of this sort is unstable in that an infinites- 
imal anisotropy is enough to trigger a run-away towards flavor equipartition. The "self-maintained 
coherence" of a perfectly isotropic gas gives way to "self-induced decoherence." 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Neutrino oscillations between two flavors in vacuum 
are governed by the frequency 
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where E is the energy of a given mode. Therefore, if 
the energy spectrum is broad, the energy dependence of 
the oscillation frequency quickly leads to kinematical de- 
coherence, i.e., along a neutrino beam the overall flavor 
content quickly approaches an average value. 

The situation changes radically when neutrinos them- 
selves provide a significant refractive effect, leading to 
collective oscillation modes [H, 0, i, H i, 1, 0, S, i, 
Tol . [Til [HI that can be of practical interest in the 
14 , IH fill o r in core-collapse super- 
Defining the 



early univers e [Vc 



novae 0, HI [llMMlfll, HSU 

parameter 



Li = V2G F n v , (2) 

the neutrino gas is "dense" when ui < [i, i.e., when a typ- 
ical neutrino-neutrino interaction energy exceeds the en- 
ergy corresponding to the vacuum oscillation frequency. 
When this condition is satisfied, collective effects are im- 
portant, even if the ordinary matter effect is much larger 
than that from the neutrino- neutrino interaction [23|, [2g] . 
One characteristic feature of collective oscillations is the 
phenomenon of "self-maintained coherence" [6| , meaning 
that all modes oscillate with the same frequency even 
though the energy spectrum may be broad. 

It was recently stressed, however, that this is not 
the complete story [ill, Hj]- Perhaps the most inter- 
esting case for collective effects is provided by neutrinos 
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streaming off a supernova core, a situation that is far 
from isotropic. The current-current nature of the weak- 
interaction Hamiltonian implies that the interaction en- 
ergy between particles of momenta p and q is propor- 
tional to (1 — v p • v q ) where v p = p/E is the velocity. 
If the medium is isotropic, the v p • v q term averages to 
zero, but if there is a net current, test particles moving 
in different directions experience different refractive ef- 
fects. Therefore, neutrinos moving in a background with 
a net current will decohere between different directions of 
motion, even if the energy spectrum is monochromatic. 

To avoid confusion about terminology we stress that 
we always mean "kinematical decoherence" when we 
say "decoherence." Different modes oscillate differently, 
leading to "de-phasing" and thus to the depolarization 
of the overall ensemble. If we use the common language 
of polarization vectors P p for each mode p, then the 
length of each P p is conserved, whereas the length of 
the overall polarization vector P = Y] P p shrinks (kine- 
matical decoherence) or is conserved (kinematical coher- 
ence). The effect of dynamical decoherence, caused by 
collisions among the neutrinos or with a thermal back- 
ground medium, is that each individual polarization vec- 
tor P p shrinks, i.e., neutrinos in individual modes cannot 
be represented as pure states. This effect, relevant for 
open quantum systems, does not occur in our case where 
oscillations are the only form of evolution. 

The multi-angle decoherence effect becomes nontrivial 
in the most interesting case when the "background cur- 
rent" consists of the neutrinos themselves as for neutri- 
nos streaming off a supernova core. Numerical examples 
reveal significant decoherence effects, but on the other 
hand they also show collective modes very similar to the 
isotropic ("single-angle") case [13] • The overall behav- 
ior is determined by a complicated interplay between the 
collective evolution and decoherence. 

Rather than trying to understand the supernova case 
in its full complexity, we here take the opposite ap- 
proach and study the simplest example that shows de- 
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coherence caused by the neutrino-neutrino multi-angle 
effect. Therefore, we consider a dense neutrino gas that 
is monochromatic, symmetric (equal v and v densities), 
and homogeneous, but not isotropic. Even this simple 
model has a rich phenomenology that helps one to de- 
velop a better understanding of the full problem. On the 
other hand, it is simple enough that the most puzzling 
aspects of its behavior are analytically tractable. 

A symmetric vv gas oscillates in the "bipolar mode" 
where pairs of neutrinos of a given flavor coherently os- 
cillate into the other flavor and back with the "bipolar 
frequency" @, % M M 



(3) 



As an example we consider a dense neutrino gas initially 
consisting of v e v e and assume that the mixing angle 
with another flavor is small, a situation that could rep- 
resent 13-oscillations. We characterize the overall flavor 
content of the ensemble in the usual way by a polarization 
vector where the positive z-direction corresponds to the 
i/ e -flavor whereas the mass direction is defined by a unit 
vector that we call B. It is tilted against the z-direction 
by 20 and we always choose it to lie in the x-z-plane 
with B z < for the normal hierarchy and B z > for the 
inverted hierarchy. 

In vacuum, the polarization vectors precess around B 
with frequency to. In a dense neutrino gas with 8 < 1, 
their motion is largely confined to the x-z-plane where 
they perform pendular motions with frequency k and 
where — B is the "force direction." The maximum excur- 
sion is 20, indicated by the dotted lines in Fig. [JJ where 
we show the trajectory in the x-z-plane. For the normal 
hierarchy, the pendulum swings between the dotted lines. 
For the inverted hierarchy (bottom), it performs almost 
full-circle oscillations as indicated by the diamonds. 

If the gas is not isotropic, the picture changes com- 
pletely. Taking the "half-isotropic" case as a generic ex- 
ample where only one half-space of momentum modes is 
occupied (neutrinos streaming off a plane surface), we 
show in Fig. [1] the trajectory of the average polarization 
vector. For both hierarchies, its length shrinks to zero 
so that the final state of the ensemble is an equal mix- 
ture of both flavors. In other words, the ensemble quickly 
decoheres, independently of the mass hierarchy. This be- 
havior is perhaps intuitive because we have assumed a 
large neutrino current. 

However, this effect is self-accelerating in that it is trig- 
gered even by the smallest nonzero anisotropy. Whereas 
for perfect isotropy the neutrino-neutrino interaction pro- 
vides for self-maintained coherence between modes of dif- 
ferent energy, even an infinitesimal anisotropy has the op- 
posite effect of causing self-induced decoherence between 
modes with different propagation directions. The aver- 
age flavor content of the ensemble always approaches an 
equal mixture of both flavors. This instability is a non- 
linear effect caused by the neutrino-neutrino interaction. 
In other words, the pendular motion of the isotropic en- 
semble is an unstable limit cycle, whereas the equal flavor 
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FIG. 1: Trajectory of the overall polarization vector in the 
"half-isotropic case" for the normal hierarchy (top) and in- 
verted hierarchy (bottom) with sin26> = 0.1. The arrow is 
the B direction. In the isotropic case (no decoherence), the 
dotted lines represent the largest pendulum excursion of 2G 
on either side of B. The diamonds in the lower panel represent 
the pendular motion in time intervals of 0.01 (2tt/k). 



mixture is a stable fixed point. 

Our goal is to explain this behavior. As a first step 
we formulate in Sect. [Til the equations of motion for the 
multi-angle system in terms of the multipole components 
of the neutrino angular distribution. The special role of 
the flux (the dipole term of the angular distribution) be- 
comes manifest. We begin with the full ./V-flavor problem 
and then consider two flavors in terms of polarization vec- 
tors. In Sect. IIIII we study the role of a large anisotropy 
in causing decoherence. Our main result is developed in 
Sect. lIVI where we show that a vv gas is unstable in flavor 
space and always decoheres if there is an infinitesimally 
small initial anisotropy. We conclude in Sec. [V] with a 
summary and outlook. 
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II. MULTI-ANGLE EQUATIONS OF MOTION 
A. General equations of motion 

A statistical ensemble of unmixed neutrinos is repre- 
sented by the occupation numbers / p = (a p et p ) for each 
momentum mode p, where ap and a p are the relevant 
creation and annihilation operators and (. . .) is the ex- 
pectation value. The corresponding expression for an- 
tineutrinos is / p = (a p a p ), where here and henceforth 
overbarred quantities refer to antiparticles. 

For several flavors, the occupation numbers are _gener- 
alized to "matrices of densities" in flavor space 30, 31. 

P )ij = (ajai) P and (g p )ij = (a] a,j) p (4) 
that really are matrices of occupation numbers. The re- 



versed order of the flavor indices i and j in the r.h.s. for 
antineutrinos ensures that under a flavor transformation, 
g p — > U g p U\ antineutrinos transform in the same way, 
g p — > UgpU^. Sums and differences of g p and g p then 
transform consistently. For the seemingly intuitive equal 
order of flavor indices that is almost universally used in 
the literature, g p , g p , gp and gp all appear in the equa- 
tions instead of g p and g p alone. 

Analogous matrices (£ p )ij can be defined for the 
charged leptons. 

We assume that the neutrino ensemble is completely 
described by these one-particle occupation number ma- 
trices, i.e^, that genuine many-body effects play no 
role [13, HH, [111- In this case flavor oscillations of an 
ensemble of neutrinos and antineutrinos are described by 
an equation of motion for each mode [3(| HH, [H, [Hj] 



id t g p = +[n p ,g p ] + V2G F 
i<9*f? P = - [ftp, Qp] + V2G F 



d 3 q 



(2tt) 3 
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Vp),£>p 



v q - v p ),g p 



(5) 



where [•, •] is a commutator and Gf is the Fermi constant. For ultrarelativistic neutrinos, the matrix of vacuum 
oscillation frequencies, expressed in the mass basis, is il p — diag(mf , m|, m§)/2p with p = |p|. Further, v p = p/E p 
is the velocity of a particle (neutrino or charged lepton) with momentum p. 

The total conserved energy of the neutrino ensemble is an important quantity for understanding its evolution [23l. |26| . 
We find for the energy density 



e = 



J J^ n p( g P + ^ + ^ GF J J2nf j (2^3 (^-^(^-^(i-Vq-Vp) 



V2G F 



(27T) 



3 (Qp ~ Qp) 



V2G F 



d 3 T 



(Qp ~ ?p)v P 



(6) 



This quantity actually represents the energy shift caused 
by the neutrino masses and by neutrino interactions. It 
is straightfoward to show that indeed e — by taking the 
time derivative on the r.h.s. of Eq. ([6|), inserting the equa- 
tions of motion Eq. (JSJ) , and using cyclic permutations of 
matrices under the trace. 

For completeness we mention that one can also define 
the entropy density fjjlj 

/d 3 p 
^—3 Tr[g p \n{g p ) + (1 - g p ) ln(l - g p ) 

+ Q P HQp) + (1 ~ §p) ln(l - Qp)] • (7) 

This expression is well defined because g p and 1 — g p are 
positive semi-definite matrices. In our case where the 
equation of motion for each g p is of the form idtg p = 
[H,g p ], the entropy density is conserved. Oscillations 



alone do not lead to a loss of information. Kincmatical 
decoherence does not lead to an increase of entropy. 



B. Axial symmetry 

Since we wish to study the simplest example that shows 
nontrivial multi-angle effects, we now restrict the neu- 
trino energy distribution to be monochromatic and the 
geometry to axial symmetry. We then integrate the ma- 
trices of densities over all variables except u = cos# p 
where the angle is relative to the direction of symmetry. 
Therefore, we consider the matrices 

1 f dpp 2 df 
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implying the normalization 



1 



+ 1 



go = -J -w?** = U dU6u ' (9) 

For convenience we have arbitrarily normalized the g u 
matrices to the neutrino density. 

We further define the matrix representing the particle 
flux along the direction of symmetry, 



Ql 



du q u v v 



(10) 



where v u is the velocity along the symmetry direction. 
For relativistic particles v u = cos 6 U = u. 

For charged leptons we proceed in the same way, except 
that we normalize to the electron density n e . 

We denote with to the matrix of vacuum oscillation fre- 
quencies for our fixed neutrino energy E. The equations 
of motion then simplify to 

ig u = [to,g u ] + (£ ~ ()o + h(q ~ o) , Qu] 

- [X(i-l) l +fi(g-g) 1 ,ug u \ ,(11) 

where A = \/2GFn e . The same equation applies for p u 
except for a sign change of the vacuum oscillation term. 

Evidently the isotropic part of the medium (index 0) 
affects all modes in the same way and is ultimately re- 
sponsible for the phenomenon of self-maintained coher- 
ence. On the other hand, the flux term (index 1) involves 
a factor u = cos# for every mode u. Even in the absence 
of a neutrino-neutrino term, a charged-lepton flux alone 
causes a trivial multi-angle decoherence effect. 

The conserved energy of the axially symmetric neu- 
trino ensemble is 

E = Tr[fi(f> + £o)] 

+ A Tr [(£ - £) (g - g) -{I- l) x {g - g),] 

+ ^[{g--g)l-(g-g)\] . (12) 

Because of our normalization of the g u matrices the quan- 
tity E = e/n v is the energy per v. 



C. Expansion in Legendre polynomials 

The structure of these equations becomes more trans- 
parent if we expand the angular dependence of the matri- 
ces in Legendre polynomials. (Had we not assumed axial 
symmetry, spherical harmonics would be the appropriate 
basis.) The first few Legendre polynomials are 



L (u) = 1 , 
Li(u) = u, 
L 2 (u) = \{Zu 2 -l). 



We thus define 



+i 



dug u L n (u) 



(13) 
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The normalisation 
-l 



duL m {u)L n (u) 



2n + l 



implies that the original function is 

oo 

](n + \)g n L n {u) . 



g u 



71=0 



(15) 



(16) 



The previously defined overall density go and the flux g\ 
are but the zeroth and first case of Eq. (fT4")) . 

Multiplying both sides of Eq. (fTTj) with L n (u) and in- 
tegrating over du leads to the equation of motion 



ig n = [to, g n ] + ju [(q - g) , Q n ] 

(Q- q)i, S^ x duuQ u L n {u) 



A 4 



(17) 



Here and henceforth we no longer show the ordinary mat- 
ter term. Its structure is similar to the neutrino-neutrino 
term so that it is easily reinstated. 

With the expansion Eq. p6[) the remaining integral is 



duug u L n (u) 



m=0 



{m + \)g m l duuL m (u)L n (u) . (18) 



We use the "dipole matrix element" of the Legendre poly- 
nomials 



+ \ T ( \T ( \ 2(m+l) 

duuL m (u)L n {u) = — — — — 

_i (2m + l)(2m + 3) 

2m 



(2m- l)(2m+ 1) 



&m + 1 , n 



With this result the equations of motion are 

ig n = [to,Q n ] + fJ,[(Q- Q) 0> Q n ] 



(19) 



(20) 



where 



2n 



2n + l 
2(n + 1) 
2n + 1 



= 1 - 
= 1 



1 



2n + 1 : 
1 

2n+ 1 ' 



(21) 



Note that the equation for go is consistent even though 
the quantity g-\ is not defined because the coefficient 
a n = n l (2« + 1) = for n = 0. 

The equations of motion Eq. (|2H)l together with the 
corresponding equations for g n (sign change for the vac- 
uum term) form a closed set of equations. If we also 
include the ordinary matter term they are equivalent to 
Eq. (jlip for the momentum-space matrices g u . 
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One important difference is that Tr(p^) is conserved, 
reflecting the absence of decoherence for individual mo- 
mentum modes, whereas Tr(f?^) is not in general con- 
served, reflecting the effect of kinematical decoherence. 
We note that in general <9tTr(A 2 ) — when the equation 
of motion is of the form i9(A = [A, H]. This is the case 
for g u , but not for the multipole matrices g n . 

Of course, in a numerical implementation the series g n 
has to be truncated at some value n max , leading to lim- 
ited angular resolution. This is analogous to the coarse 
graining of phase space required for the p u where one 
needs to use discrete angular bins of nonzero width. 

The equations of motion for the flux terms are special 
in that they involve one power of Q\ or g~\ in each term 
of the equation. Therefore, if initially there is no flux 
term (gi = g± = 0), none will develop. In this case the 
equations for g$ and g~Q form a closed set, describing the 
dynamics of the "flavor pendulum" studied in Ref . (2(| . 
In addition, the higher multipoles g n with n > 2, if ini- 
tially nonzero, will simply oscillate under the action of 
the vacuum term and of the density term (g — g)o- 



D. Diffusion equation in multipole space 

The interpretation of the equations of motion in mul- 
tipole space becomes clearer if we observe that for some 
function f(x) one has to second order 

fix + Ax) = f(x) + f'(x)Ax + \ fix) Ax 2 , 

f(x - Ax) = f(x) - f'(x)Ax + \ f"{x)Ax 2 . (22) 

Taking the sum and difference of these equations provides 



f(x + Ax)+f(x- Ax) 
f(x + Ax)-f(x-Ax) 

With Eq. (f2"T) we may write 
flngn-i +b n g n +i — 2g r , 



2n+ 1 



2f(x) + f"(x)Ax 2 , 
2f(x)Ax. (23) 



g' n An + g^An 2 , (24) 



where we interpret a prime as a derivative with respect 
to n that is now viewed as a continuous variable. If we 
further observe that in the discrete equation An = 1 wc 
may write Eq. (|20[) in the continuous form 



\g n = [Cl,g n ] +n[(g- g) - {g- g)x,g n ] 



(Q- q)i, 



2n + 1 



,(25) 



Therefore, the lowest multipole (monopole) is equivalent 
to the pendulum in flavor space [261 ] , the high multipoles 
are akin to a continuous medium carrying flavor waves. 
Equation (|25[) is a diffusion equation. The periodic ex- 
citation caused by the flavor pendulum then diffuses to 
higher multipoles (smaller scales) . The flavor waves as a 
function of neutrino direction visible in numerical simu- 
lations [25j are the result of this process. 



E. Particle vs. lepton number 

In the equations of motion the difference of the parti- 
cle and antiparticle densities, g p — g p , plays a special role 
in that it is this matrix of net lepton number densities 
that plays the role of a potential for the other neutrino 
modes. As noted by several authors, it is sometimes use- 
ful to write the equations of motion in terms of the sum 
and difference of g p and g p instead of these matrices 
themselves. We therefore define the matrix of particle 
densities (sum S) and the matrix of lepton number den- 
sities (difference D) by virtue of 



D r 



Qp ~ Qp 



(26) 



Analogous definitions apply to the case of axial symme- 
try, S u and D„, and their multipole expansion S„ and D n . 

The equations of motion for the particle and lepton 
number matrices are found by adding and subtracting 
the two lines of Eq. ([5]) and the equivalent equations for 
axial symmetry or the multipoles. We find explicitly for 
the multipoles 



iS„ = [CI, D„] + ju[D ,S„] - n 
iD n = [Cl,S n ] +^[D ,D„] - [i 



Di 



Di, 



anD n _i + 6„D n+1 



(27) 



One can easily restore the ordinary matter term because 
it has the same structure as the neutrino-neutrino term 
with the substitution u — > A and Dc^i — > Lo.i, the latter 
being the charged-lepton matrix of net lepton-numbcr 
densities and the corresponding flux, respectively. 

The commutator structure of the r.h.s. of Eq. (|27p im- 
plies that its trace vanishes. Therefore, the particle and 
lepton numbers for all multipoles, Tr(S„) and Tr(D„), 
are conserved. 

With the expressions Eq. (f2"Tj) for a n and b n the lowest- 
order multipole equation is 

iS = [0,D ]+/i[D ,So]-M[Di,Si], 

iD = [n,S ]. (28) 

One immediate consequence of the second equation is 

d t Ti(ClD ) = , (29) 

i.e., in the mass basis a certain combination of flavor- 
lcpton numbers is conserved. 

The explicit equation for the flux term is 



iSj = [fi,D 1 ]+^[D ,S 1 ] + |[S + 2S 2 ,D 1 ] 
iDx = [n,Si] + ^[D -D a ,D 1 ]. 



(30) 



As stated earlier, these equations are special in that they 
involve one power of Si or Di in each term. Therefore, 
Si and Di vanish at all times if they vanish initially. 
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In this case, the monopole terms (overall densities) 
form a closed set of equations. Defining 



Q 
D 



So- 
Do, 



/' 



the equations of motion take on the simple form 

iQ = /i[D,Q], 
iD = [«,Q]. 



(31) 



(32) 



Conserved quantities are Tr(Q 2 ), Tr(DQ), Tr(DQ), and 
Tr(f2Q + ptD 2 /2). Up to a constant, the last quantity is 
the energy per neutrino. For two flavors, these simpli- 
fications are at the origin of the pendulum analogy [261 ] 
where Tr(Q 2 ) corresponds to its length, Tr(DQ) its spin, 
and Tr(Df2) to the orbital angular momentum around the 
force direction (the mass direction in flavor space). 

Returning to the general case with arbitrary initial 
conditions, the energy per neutrino is 



E = Tr 



ttSo + f (D 2 -D?) 



(33) 



Energy conservation dtE = is easily verified by using 
the equations of motion for So, Dq and Di. 



Two-flavor case 



is represented by 



d t Tr(A) =0 and 9 ( A = BxC. 



(36) 



Moreover, 



Tr(AB) 



Tr(A) Tr(B) + A B 



(37) 



for any two matrices A and B. 
As usual, we denote with P r 



and P p the polariza- 



tion vectors representing g p and g p , respectively, and 
analogous expressions for g u — > P u and g n — » P„. The 
particle and lepton-number matrices are represented by 
S — * S and D — * D, in each case with the subscripts p, u 
or n for the different variables. Finally, we represent the 
matrix of oscillation frequencies in the form 



n 



Tr(f2) + wB • tr 



(38) 



where 



Tr(O) 



2E 



and u> 



2E 



(39) 



In the two-flavor case where the g p are Hermitian 2x2 
matrices, a representation in terms of polarization vectors 
is useful. A Hermitian 2x2 matrix A is represented as 



A = 



Tr(A) + A • tr 



(34) 



where tr is the vector of Pauli matrices and A the polar- 
ization vector. 

The commutation relations of the Pauli matrices imply 
that an equation of motion of the form 



iftA = [B, C] 



(35) 



We have defined the oscillation frequency u> as a positive 
number. The "magnetic field" is a unit vector that in 
the interaction basis is explicitly B = (sin 29, 0, cos 20) 
where is the vacuum mixing angle. We have arbitrar- 
ily chosen its y-component to vanish, corresponding to 
a choice of overall phase of the mixing matrix. In the 
following, when we say that "the mixing angle is small," 
we mean that | sin 20 1 <C 1. The sign of B x is physically 
irrelevant whereas B z < represents the normal, B z > 
the inverted hierarchy. 

The equations of motion for the particle- and lepton- 
number polarization vectors are explicitly 



S„ = ljB x D„ + fiD xS„-^DiX (a„S„_i + b n S n+1 ) , 

D„ = ujB x S n + ^D x D„ - | D! x (a n D n _i + b n T> n+1 ) . (40) 



I 

The two lowest-order equations are explicitly and 

Si = uBxDi + fiBo x Si + | fj, (S + 2S 2 ) x Di , 
S = uBxD + Li&o x S - /i Di x Si , Di = cjB x Si + § fi (D - D 2 ) x Di . (42) 

D = ujBx S (41) 
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The conserved energy Eq. p3|) becomes, with the help 
of Eq. (J37J), 

E = i{Tr(ft)Tr(S )+^B-S 

+ | [Tr(Do) 2 - Tr (D^ 2 + D 2 - D 2 ] } . (43) 

We have already noted that the individual traces Tr(So), 
Tr(D ) and Tr(Di) are conserved so that 1 



E =2 



;B • So + | (D 2 - D 2 ) 



(44) 



is the conserved energy per neutrino, up to an irrelevant 
constant. We use the notation 



E u 



> 



o 



B-S ) 



(45) 



for what is the "potential energy" in the pendulum anal- 
ogy. Its minimum is chosen at zero. The "kinetic energy" 
(neutrino-neutrino interaction energy) is denoted by 



(46) 



where we have introduced 

E Q = 



^D 2 



! 



for the monopole and dipole contributions. 



(47) 



III. ROLE OF ANISOTROPY 

A. Isotropic case 

In order to study the impact of kinematical decoher- 
ence in the simplest nontrivial case we study a homoge- 
neous ensemble consisting of equal numbers of neutrinos 
and antineutrinos that are all initially in the electron fla- 
vor. Before turning to issues of decoherence, we briefly 
recall the evolution of the isotropic case. To this end we 
show in Fig. [2] in the upper panels the z-component of the 
overall polarization vector So = Po + Po for the inverted 
and normal hierarchy. The v e or v e survival probability is 



prob(^ e — * v e ) = prob(^ e — > v e ) 



So 2 + 2 



(48) 



because in our normalization |So| =2. In the lower pan- 
els we show the corresponding energy components. 

For the normal hierarchy (left) the motion is a small- 
excursion harmonic oscillation with the bipolar fre- 
quency k. We show one full cycle of oscillation. The 



1 The overall factor 1/2 is missing in the equivalent expression in 
our previous paper |2(| . 



potential energy E^ (solid line in the lower panel) is mea- 
sured against the mass direction so that it is back at its 
maximum after a half-period as behooves an ordinary 
pendulum. The dashed line shows the "kinetic pendu- 
lum energy" E^. The isotropy of the ensemble implies 
E\ = and thus — E . 

For the inverted hierarchy (right) the motion is that 
of an inverted oscillator. We show approximately one 
full cycle that for the chosen mixing angle lasts approx- 
imately 3 times the bipolar period 2tt/k. The polariza- 
tion vector starts in the positive z direction, swings al- 
most one full circle and almost arrives back at the verti- 
cal direction (maximum in the middle of the plot), then 
swings back, reaching its original vertical position, and 
then starts over again. The maximum potential energy is 
2w because the energy is normalized to the total energy 
per v and because there are equal numbers of neutrinos 
and antineutrinos, i.e., two particles per v. 

Reducing the mixing angle has the effect of reduc- 
ing the oscillation amplitude for the normal hierarchy, 
whereas it increases logarithmically the duration of the 
"plateau phases" for the inverted hierarchy. For quan- 
titative discussions of the isotropic case in terms of the 
pendulum analogy see Ref. [26| . 



B. Half-isotropic case 

Next we consider the same example, but assume a large 
degree of anisotropy where only one half-space of mo- 
mentum modes is occupied. In Fig. [3] we show, in the 
upper panels, the evolution of 5q. Moreover, we show 
the length of the overall polarization vector, |So|, as a 
thin dashed line. Both in the normal hierarchy (left) and 
the inverted hierarchy (right), an equal mixture of the 
two flavors is quickly achieved. In the isotropic case the 
length of Po and Po is conserved and that of So is approx- 
imately conserved, up to corrections of order uj/fi <C 1. 
Here, the lengths of Po, Po and So shrink to zero, reflect- 
ing kinematical decoherence. Of course, Sq can become 
small or zero without its length shrinking as during the 
first swing in the inverted hierarchy case. It is the length 
of the polarization vector, not its z-component, that is a 
measure of decoherence. Of course, in the normal hierar- 
chy, So performs only small-excursion oscillations so that 
a significant change of its z-component can be achieved 
only by a change of its length. In this case |So| and Sq 
are almost identical. 

The evolution of So was also illustrated in Fig. [T] where 
we showed its trajectory in the x-z-plane. For our cho- 
sen geometry where B is in the x-z-plane, Sq = at all 
times. For the inverted hierarchy (bottom) we indicate 
the isotropic-gas trajectory with diamonds at time inter- 
vals of 0.01 (2tt/k) where the motion starts in the vertical 
position. In the half-isotropic case, the particlc-numbcr 
polarization vector So spirals in to a position close to the 
origin. The final offset is very small and depends on the 
magnitude of for which we have chosen 10~ 5 . In 



8 




FIG. 2: Evolution of a homogeneous and isotropic ensemble of equal numbers of neutrinos and antineutrinos, initially in the 
electron flavor, with a vacuum mixing angle sin 20 = 0.1. Top: Polarization vector S§- Bottom: Energy components (solid) 
and Efj, (dashed). Left: Normal hierarchy. Right: Inverted hierarchy. 




FIG. 3: Same as Fig. [2] now for a half-isotropic neutrino gas. Top: Polarization vector So (solid) and the length | So | (dashed). 
Bottom: Energy components (solid), _B M = Eq + E\ (dashed), Eo (thin dotted) and E\ (thin dot-dashed). 
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the normal hierarchy (top), the isotropic-case motion is 
a small-excursion harmonic oscillation between the dot- 
ted lines. In the half-isotropic case, the evolution begins 
from the vertical position. A full oscillation back to the 
vertical position is performed before the shrinking of So 
becomes noticeable. 

In the lower panels of Fig. [3] we show the evolution 
of the different energy components. The simpler case is 
the inverted hierarchy (right) where the "potential" and 
"kinetic" energies begin to oscillate as in the isotropic 
case of Fig. [5J On the time scale of a few bipolar oscil- 
lation periods, the two components essentially equiparti- 
tion, although asymptotically a small offset remains that 
depends on u / fi. Moreover, the neutrino-neutrino energy 
( "kinetic energy" ) now develops a nonvanishing flux term 
E\ = —/.iDl/4 that inevitably is negative. 

The normal hierarchy (left) is initially similar in that 
and E^ oscillate as for a pendulum, even though this 
motion is not visible on the scale of the plot. As deco- 
herence sets in, a qualitatively different mode of behavior 
obtains in that the neutrino-neutrino energy £?„ is dom- 
inated almost entirely by the negative E\ whereas Eq 
now is subdominant. All individual polarization vectors 
of all modes start aligned with the z-direction, i.e., al- 
most aligned with the force direction since the mixing 
angle is small. Therefore, the initial E u is near its mini- 
mum. If the overall polarization vector So is supposed to 
shrink, the individual polarization vectors must develop 
significant deviations from each other and thus the po- 
tential energy must increase. Energy conservation then 
dictates that E^ becomes negative. 

One important conclusion is that the angular depen- 
dence of the neutrino-neutrino interaction alone is not 
enough to cause kinematical decoherence, but its abso- 
lute sign is also crucial. If E\ were not negative, energy 
conservation would prevent significant decoherence for 
the normal hierarchy. Changing this sign in a numerical 
example indeed reveals the absence of decoherence for the 
normal hierarchy, but no significant change of behavior 
for the inverted case. In the real world there is no choice 
about this sign. It derives from the negative sign of the 
spatial part in the neutrino current-current interaction, 
i.e., it is the negative sign inherent in the Lorentz metric. 
From the equations of motion in the form Eq. (|40[) one 
would have never guessed that the absolute sign of the 
term proportional to Di plays a crucial role. 



C. Evolution of multipoles 

It is also instructive to study the evolution of the higher 
multipole components S„ . We consider the same example 
as in the previous section. For a half-isotropic gas the 
initial relative length of the multipole components is 



S 



Si 
So 

S2m 

So 



1 

2 " 
0, 



2m+l 

~So~ 



1 ^ ( 1 

; k=l v 



2(m+ 1) 
(-l) r 

2(m+l) r(|)r(m + i) ' 



(49) 



where S n = |S n | and all of them are initially oriented 
in the ^-direction. Here, m is an integer, i.e., the even 
multipoles vanish initially except for So- 

Instead of using these initial values, we calculate the 
evolution for a case where all multipoles vanish initially 
except for So = 2 and S\ = 1, i.e., we use the same 
density and flux terms as in the half-isotropic case. The 
evolution is very similar. We show the evolution of the 
first 21 multipoles in Fig. [4] for the inverted hierarchy; 
the picture is similar for the normal hierarchy. The ini- 
tial horizontal part of the curves corresponds to their 
zero initial value so that the offset by 0.1 vertical units 
between the curves is directly apparent. 
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FIG. 4: Evolution of |S„| for the first 21 multipoles in 
the quasi-halfisotropic example with inverted hierarchy. The 
curves are offset from each other by 0.1 vertical units. 
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FIG. 5: |S n | (solid) and S% (dotted) as a function of n at 
t — 12 for the example of Fig. [3] 
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Initially, So and Si are large, but quickly decay away 
by decohcrcnce, whereas the higher multipoles vanish ini- 
tially and get excited one after another, but then decrease 
again. The length of each multipole sports "wiggles," 
i.e., the evolution is not monotonic. Moreover, they have 
large relative angles in the x-z plane (not shown here), 
i.e., the spread of the inital excitation to larger and larger 
multipoles is far from simple in detail, but the overall 
process is as expected. In Fig. [5] we show |S„| and S* 
as a function of n at t = 12 for the example of Fig. 2] 
At t = 12, multipoles n > 60 are not yet excited. The 
nearly linear increase of the phases between individual 
polarization vectors corresponds here to a diffusion of 
the "multipole wave" to larger n with approximately con- 
stant speed, n ~ nt. 

Qualitatively, this can be understood as follows. The 
second line in Eq. (HP)) together with the analogue of 
Eq. (|24|) implies that 



D„ ~ wB x S„ — ^Di x 



d: 



2n 

/i(Do-Di) 



f 1 
x D, 



D" 

2 



(50) 



The last term does not contribute to the growth of 
|D„| and can be ignored. Approximate equipartition 
between potential and kinetic energies in Eq. (|44p im- 
plies |D | ~ |Di| ~ 2(w/ ( u) 1 / 2 , as indeed was observed 
in Fig. [3] The second term in Eq. ([50)) then repre- 
sents a combination of a drift in n— space with veloc- 
ity ~ K;/(v2n) and diffusion with diffusion coefficient 
~ k/\/2. Up to a factor < -1 / 2 this gives rise to an expo- 
nential factor ~ exp[— n /(2v / 2k*)]- On the other hand, 
since |S„| ~ (fj,/u>) ' 2 \D n \ for n > 2, the first term in 
Eq. (|50p can give rise to exponential growth whenever 
(B x S n ) • D„ > 0. This growth will be non-monotonic, 
but can roughly be estimated as cx exp(nt/2) before D„ 
saturates. Using perturbation theory, we will find a simi- 
lar growth rate for Di at early times in Sect. lIVBl Com- 
bining these two factors, one sees that at time t the wave 
front of D„, and thus also of S n , should be located at 
n ~ Kt, which, according to Figs. |4] and [5] indeed it is 
within about 20%. The slow decrease of S„ and D„ at 
late times can be interpreted as due to the t -1 / 2 factor 
in the diffusive behavior. 



D. A simple model of decoherence 
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FIG. 6: Evolution of P§ (solid) and |P ] (dashed) of our toy 
model with Au — 1. 
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FIG. 7: Trajectory of Po in the 3>z-plane for our toy model 
of decoherence with An = 1. 



If the frequencies are spread over 1 < u < 1 + Aw, then 
the overall polarization vector evolves as 



pz 



p.r 



1+A ", . , sin[(l + Au) t] - sin(t) 

i t 

1+Am , . . . cos(t) - cos[(l + Am) t] 
du sm(ut) — , 

[(^) 2 + (P -) 2 ] 1/2 ^ 2|Sin(A ^ /2)l . (52) 



The overall behavior of the inverted-hierarchy case is 
qualitatively easy to understand. To this end we consider 
an ensemble of polarization vectors P M , initially oriented 
in the z-direction, that rotate around the y-axis, each 
with a frequency u, so that 



pz 
px 



cos(ut) , 
sin(irf) . 



(51) 



Both u and t are normalized to be dimensionless. 



For Au = 1 we show in Fig. [6] the evolution of P§ (solid 
line) and the length |Pq| (dashed line). In Fig.[7]we show 
the trajectory of Po in the x-z-plane. 

Unsurprisingly, this simple model represents qualita- 
tively the features of the inverted-hierarchy case. In par- 
ticular, the trajectory in the x-z-plane is not a true spiral, 
but there are many crossings of the origin as the length 
of the overall polarization vector shrinks. Its length does 
not shrink monotonically. Its envelope decreases as a 
power law t , not exponentially. 
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This example also shows that it is not trivial to de- 
fine a useful measure of kinematical decoherence. The 
z-component of the polarization vector is not useful be- 
cause Po can tilt so that its z-component vanishes, i.e., 
both flavors are present with the same probability, yet 
the flavor content of the ensemble is perfectly coherent. 
The length of Po is a much better measure in analogy to 
dynamical decoherence of a single mode where the length 
of P p is a measure of dynamical decoherence; unit length 
would correspond to a pure state. In our case, the length 
of the individual P p is conserved, whereas the length of 
the total Po shrinks, but not monotonically. 

It is not evident if there is an "entropy" measure that 
evolves monotonically as a consequence of kinematical 
decoherence. We note, however, that the origin of kine- 
matical decoherence is the differential evolution between 
neighboring polarization vectors. In our toy example, the 
angle between neighboring polarization vectors grows lin- 
early as Am t and thus is a measure of decoherence. 

Motivated by this observation we can define the "wind- 
ing number" 



N = — 
2n 



- 1 \dP u /du\ 

UU — — ; 



(53) 



In the symmetric case where the evolution of the polar- 
ization vectors is essentially restricted to the s-z-planc, 
this quantity tells us the number of windings around the 
y-direction of the full ensemble of polarization vectors. In 
Fig. [5] we show the evolution of this quantity for the nor- 
mal and inverted hierarchy examples of Fig. [3] N counts 
how often the length of So shrinks to zero. 
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FIG. 8: Evolution of the "winding number" as defined in 
Eq. (|53p for the half-isotropic case of Fig. [3] Normal hierarchy 
(solid) and inverted hierarchy (dashed). 



is represented by a coarse-grained (binned) ensemble of 
polarization vectors or, in multipole space, when the se- 
ries is truncated at some multipole n max . Binning or 
truncation are unavoidable in numerical treatments and 
represent equivalent approximations. 2 

If the phase differences evolve linearly between differ- 
ent polarization vectors as in the toy example of the pre- 
vious section, then the initial state will recur when each 
polarization vector has turned around the y-axis by 2tt 
relative to its neighbors. Put another way, if we use n po i 
bins to represent the polarization vectors, we expect re- 
currence to begin when N = (n po \ — l)/2. We demon- 
strate this effect in Fig. |H] for the same example as in 
Fig. [3l normal hierarchy, using n po \ = 51. Indeed, the 
winding number increases until approximately N = 25 
and then decreases back to almost zero, and so forth. 
This system is almost periodic on the recurrence time 
scale, reflecting that nonlinear effects play a subleading 
role here. Since a large flux term is present from the start, 
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FIG. 9: Evolution of |Po| and the winding number for the 
normal-hierarchy case of Fig. [3]with 51 polarization vectors. 



E. Recurrence effects 

Closely related to kinematical decoherence are recur- 
rence effects that arise when one uses a limited num- 
ber of polarization vectors, i.e., when the full ensemble 



2 In Ref. [ill a multipole expansion ("partial wave expansion") 
was proposed and dismissed with the argument that the neces- 
sary truncation would introduce uncontrolled numerical effects 
whereas coarse bins in angle space would provide limited resolu- 
tion, yet a correct solution of the nonlinear equations of motion. 
We disagree with this assessment. The same unphysical recur- 
rence effects occur both if the number of angular bins or the 
multipole order of truncation is too small. 
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kinematical decoherence is very similar to the linear toy 
example of the previous section where recurrence would 
be exact and the system would be periodic with a fre- 
quency proportional to (n po i — 1) _1 times the individual 
polarization vector's oscillation frequency. 

Using the multipoles as the primary variables, the same 
recurrence effects occur due to the truncation n max of the 
multipole series. The "multipole wave" of Fig. [5] can only 
propagate up to n max , is then reflected, and eventually 
returns back to n = 0. This is the moment of maximum 
recurrence of the overall polarization vector. 

For a more complicated example such as a realistic sim- 
ulation of the supernova case, recurrence effects may not 
be so obvious as here. Therefore, in a numerical treat- 
ment one must make sure that the largest relative angle 
developed between any two neighboring polarization vec- 
tors never grows to order unity, i.e., one should monitor 
the largest relative angle between neighboring polariza- 
tion vectors as a measure of accuracy of the calculation. 
In other words, kinematical decoherence is a property of 
the entire ensemble, but can also be a a property of indi- 
vidual bins if they are too coarse. One must make sure 
that within all individual bins kinematical decoherence 
remains negligible. 



IV. SMALL INITIAL ANISOTROPY 
A. Numerical example 

We next turn to a symmetric system that is perfectly 
isotropic except for an arbitrarily small but nonzero ini- 
tial anisotropy. We work in the multipole picture and 
assume that initially D„ = for all n, S„ = for n > 2, 
Sq — 2 and Sf = £ <C 1. Note that £ = 1 corresponds to 
a flux term equivalent to the half-isotropic case. 

As a first example we show in Fig. [10] the evolution of 
the overall polarization vector So (top) for the normal 
hierarchy (left) and the inverted hierarchy (right). The 
initial anisotropy is £ = 10~ 4 and sin 29 = 0.1. For both 
hierarchies, kinematical decoherence eventually obtains. 

In the bottom panels we show the evolution of —E± = 
(/Lt/4) ~D\ (solid line) as a proxy for Di which is the pri- 
mary quantity that causes decoherence. For the symmet- 
ric neutrino gas only the component D\ is nonzero, but 
it can be negative and can change sign so that it is diffi- 
cult to display on a logarithmic scale. We also show the 
evolution of the "potential energy" E u (dotted line). 

For the normal hierarchy, D\ grows exponentially after 
a short initial transient. While D\ is small, the potential 
energy E u performs the usual oscillations of a harmonic 
oscillator, complemented by the opposite oscillations of 
the kinetic energy E^ = Eq — E% (not shown) that is 
entirely dominated by Eq. As D\ grows, it eventually 
dominates the kinetic energy and the end state is approx- 
imately E w — +lo and E^ = —lu, with E\ dominating E^, 
In other words, the end state is exactly as in the previ- 
ous case with a large initial anisotropy. The crucial novel 



feature is that the flux term grows exponentially from a 
small value. This is a purely nonlinear effect caused by 
neutrino- neutrino interactions. For smaller mixing an- 
gles and/or other values of £, the behavior is analogous. 
The rising part of the E\ curve is the same, except that 
it shifts vertically in direct proportion to £ and sin 20. 

For the inverted hierarchy, the evolution is far more 
complicated. While D\ grows in an average sense, this 
evolution is overlaid with oscillations and in fact D\ 
changes sign at each spike in the Ei curve. Moreover, the 
evolution of the envelope of the E\ curve is not mono- 
tonic. It can perform complicated motions, with growing, 
declining, and nearly flat phases, until finally decoherence 
obtains. 

From these numerical observations we tentatively con- 
clude that an isotropic neutrino ensemble of equal densi- 
ties of neutrinos and antineutrinos is not stable. A small 
flux term triggers a run-away evolution towards kinemat- 
ical decoherence. 



B. Analytic treatment 

1. Simplified equations of motion 

As a starting point for an analytic understanding of 
this instability we use the equations of motion Eq. ([40]) 
for the vectors S„ and D„. Initially only So is of order 
one, | S 1 1 = £ -C 1, and all others vanish. As long as Di is 
sufficiently small, none of the higher multipoles can grow 
large so that we can limit our attention to the equations 
for n = and n = 1, i.e., to Eqs. (gU) and ([32]). 

We further observe that the n — equations are cou- 
pled to the higher multipoles only by the term — /1D1 x Si, 
consisting of a product of two small quantities, i.e., it is 
at least of order £ 2 . Therefore, the evolution of the near- 
isotropic ensemble is identical with that of the isotropic 
case until Di has grown sufficiently large. In other words, 
up to second order in £ the n — equation is that of the 
isotropic case. 

Therefore, what remains to be solved is the n = 1 
equation. After neglecting terms involving n > 2 all we 
need to study is 

Si = (uB + | S ) x Di + /iD x Si , 
Di = ujB x Si + ^ D x Di , (54) 

o 

where Sa(t) and Do(t) are the solutions of the unper- 
turbed pendulum equations. 

The equations simplify further in our case of equal v 
and v densities where symmetry dictates that all polar- 
ization vectors P p and P p evolve as each other's mirror 
images relative to the z-z-plane, the latter being singled 
out by our choice that B lies in that plane. In this case 
all D„ vectors are parallel to the y-axis, whereas all S„ 
vectors are confined to the x-z-plane. 
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FIG. 10: Evolution of a homogeneous and near-isotropic ensemble of equal v and v densities with sin 20 = 0.1 as in 

Figs.[2]andH 

The initial anisotropy is £ = 10~ 4 . Top: Polarization vector S§ (solid) and the length | So | (dotted). Bottom: Neutrino- neutrino 
flux energy —E\ (solid) and potential energy E u (dotted). Left: Normal hierarchy. Right: Inverted hierarchy. 



As a consequence, the Do x Di term drops out. The 
second equation then implies Di = wB x Si, and with 
the first equation yields 



Di = uiB x 



wB + | S ) x Di 



uj^iB x (D x Si) . 

(55) 

Expanding the triple product, observing that in our case 
B • D„ = 0, and using B 2 = 1 this is 



LOfJ, 



B S Di+^(B-Si)D . (56) 



Di 



We repeat that all D„ are here parallel to each other and 
to the y-axis. In the absence of neutrino-neutrino interac- 
tions (jx = 0), Di obeys a harmonic-oscillator equation. 

In a dense neutrino gas, fi > w, we may ignore the 
term proportional to uj 2 . Since during the growth phase, 
|Di| < |D | and |Si| < |S |, the last term in Eq. JSHD 
can also be neglected. A possible exponential growth is 
then accounted for by the term 



Di 



Di 



(57) 



3 2 

where we have used the bipolar oscillation frequency k 2 — 
2w/z and we note that initially |S | — 2 and that its length 
remains nearly constant if |sin20| <C 1 unless there is 
kinematical decoherence. 



2. Normal hierarchy 

If the mixing angle is small, the "normal hierarchy" is 
defined by neutrinos being essentially in the lower mass 
eigenstate or B ■ So/2 = — 1 up to corrections of order 
sin 2 29. Therefore, the critical part of the equation of 
motion for the flux term is approximately 



Di=+yDx. (58) 
The solution includes growing modes of the form 

Di cx exp ^+-^rj , (59) 

where r = nt/27t is the dimensionless time variable used 
in all of our plots, i.e., time in units of the bipolar os- 
cillation period. This behavior is numerically verified in 
the lower-left panel of Fig.QJJif we recall that the plotted 
quantity is — E± cx D 2 cx exp[(47r/vo) r]. 

Once more we note that the crucial absolute sign in 
Eq. (f58|) is traced back to the negative sign in the term 



(1 — v q ■ v p ) in Eq. ([5]). Therefore, the instability of the 
near-isotropic neutrino gas, for the normal hierarchy, is 
not just caused by different angular modes experiencing 



14 



different neutrino-neutrino effects. This would also be 
the case for a hypothetical term (1 + v q • v p ) that naively 
looks very similar, but causes a completely different over- 
all behavior. In more complicated numerical situations, 
e.g., for a full supernova simulation, it is easy to change 
this sign in order to diagnose the relevance of the term 
Eq. (|F7|) for the overall behavior of the system. 

Initially T>i ~ so that the exponential growth term 
alone is not enough. The initial evolution is dominated 
by the second line of Eq. that is initially 



x V2 



V2 



£ sin 29. 



(60) 



where D\ = |Di| and we have used |B| = 1 and Si(0) = 
(0,0, £). Here the absolute sign is not critical because it 
does not matter if Di grows in the positive or negative y- 
direction. The first factor (u/ /i) 1 / 2 sets the overall scale 
for Di. From the expression for the energy we note that 
S and D appear in the combinations w|So| and /iD§ or 
/iD 2 . The quantities S„ are of order unity, whereas the 
natural scale for D„ is (uj/fi) 1 ^ 2 . The initial evolution is 



3. Inverted hierarchy 

For the inverted hierarchy, the initial growth of D\ 
is the same (dashed line in Fig. ITT)) . After that, how- 
ever, the situation is entirely different because initially 
B • So/2 = +1. This corresponds to the pendulum start- 
ing in a nearly upright position from where it starts al- 
most full-circle oscillations. For our usual assumption 
| sin 20 1 <C 1 the term B • So is almost identical with 
Sq shown in the upper-right panel of Fig. Q7J This term 
changes sign during every swing of the pendulum. Dur- 
ing the phases when B • So > 0, Di oscillates whereas 
during the phases B • So < exponential growth obtains. 
These are the relatively short phases when the pendulum 
is oriented downward, whereas it stays upright for long 
periods if the mixing angle is small. 

We illustrate this behavior in Fig. [T2] where we show 
the evolution of D\ overlaid with that of S§. We have 
chosen an extremely small mixing angle, sin 20 = 10~ 7 , 
to obtain long plateaus for S§ where the pendulum stands 
almost still in an almost upright position. During these 
phases we have B ■ So/2 = +1 so that Eq. (|57|) implies 
harmonic oscillations of D\ with the frequency k/vo as 
borne out by the example in Fig. [T2] 



/a 1/2 

■jj) £sin29 



2tt 

71 



(61) 



where, again, r = (k/2tt) t is our usual measure of time. 
This linear growth is shown in Fig. [TT] as a dotted line 
and agrees with the numerical examples. 

Besides a transient caused by the other terms in 
the equation, the exponential growth takes over after 
a time of order 2it/k. We conclude that the initial 
value of (/Lt/a>) 1/2 Z)i for its exponential growth is of order 
£ sin (20) so that it grows to order unity within a time 
scale of order —(2it/k) ln(£sin29). 
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FIG. 11: Initial evolution of Di as discussed in the text for the 
normal hierarchy (solid) and the inverted hierarchy (dashed). 
The dotted line is the linear relationship of Eq. (|6ip . 



FIG. 12: Initial evolution of Di (solid line, arbitrary units) 
and Sq (dotted) for the inverted hierarchy with £ = 10 -4 and 
sin 29 = 10" 7 . 

The frequency is not matched to the pendulum's 

oscillation period because the duration of the "plateaus" 
depend logarithmically on sin 29. Therefore, the short 
exponential growth phases when Sq < occur at er- 
ratic instances relative to the harmonic T>i oscillation. 
This interplay explains the erratic behavior of the Di 
evolution that is apparent in Fig. [TTJ] and that obtains 
in all numerical examples. This interplay also explains 
why small changes of parameters such as sin 29 can com- 
pletely change the overall Di evolution. 

Therefore, we cannot prove if there is some specific 
combination of parameters where kinematical decoher- 
ence will not occur, although this would have to be iso- 
lated parameter points that presumably have measure 
zero in parameter space if they exist at all. 
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Except for this caveat we conclude that both for the 
normal and inverted hierarchy kinematical decoherence 
is an unavoidable consequence of the nonlinear neutrino- 
neutrino interaction terms. An infinitcsimally small, but 
nonzero, deviation from isotropy is enough to trigger an 
exponential evolution towards flavor equilibrium. There- 
fore, the pendulum solution that describes the behavior 
of a perfectly isotropic gas is merely an unstable limit 
cycle of this nonlinear system. 



V. CONCLUSIONS AND OUTLOOK 

We have investigated multi-angle kinematical decoher- 
ence effects in a dense neutrino gas consisting of equal v 
and v densities. The current-current nature of the weak- 
interaction Hamiltonian implies that a current of the 
background medium causes kinematical decoherence be- 
tween neutrinos propagating in different directions. This 
simple effect becomes entirely nontrivial in the most in- 
teresting case when the "background current" is caused 
by the neutrinos themselves. 

We have shown that a neutrino gas of this sort is not 
stable in flavor space. If one prepares the ensemble in a 
given flavor state, even a small deviation from isotropy is 
enough to trigger an exponential evolution towards flavor 
equipartition, even when the mixing angle is small. Up 
to logarithmic corrections, the time scale is determined 
by the bipolar oscillation period. We thus have to do 
with a self-induced macroscopic pair conversion process 
that proceeds much faster than ordinary pair processes 
that occur at order G|. 

With hindsight this is probably the same "speed-up 
effect" discussed in Ref. [ll(. However, the fast speed 
k = ^/2ujfi is not introduced by the multi-angle effect 
as is perhaps suggested in Ref. [ll|. Rather, k is the 
bipolar oscillation frequency that is inherent even in the 
isotropic system. For the normal hierarchy "nothing" 
seems to happen in the isotropic case, although this 
would be a misperception because the system moves with 
the speed k, but the amplitude is small if sin 20 <C 1. In 
any event, we agree with the importance of multi- angle 
effects and with the end result conjectured in Ref. [11| 
that complete flavor equipartition obtains on a time 
scale k . However, we have found no indication that 
an even faster time scale plays any role, in contrast 
to what was conjectured in Ref. [11] . 

Of course, flavor equipartition obtains only in a macro- 
scopic sense with some degree of coarse graining in phase 
space. It is the nature of kinematical decoherence that 
information is not lost and the entropy does not increase. 
On the other hand, kinematical and dynamical decoher- 
ence are not always operationally distinguishable. It is 
not necessarily possible to distinguish between a neutrino 
gas in true chemical equilibrium and one where neighbor- 
ing modes are simply "de-phased" relative to each other. 

Either way, the "self-maintained coherence" of a dense 
neutrino gas represents only an unstable limiting form 



of behavior of the perfectly isotropic case. Since no gas 
can be perfectly isotropic, self-induced decoherence is the 
generic behavior of a symmetric vv gas. 

We have only analyzed the two- flavor case. It remains 
to be studied how our results carry over to a genuine 
three-flavor situation. 

Contrary to a naive expectation, the fact that neutri- 
nos propagating in different directions experience differ- 
ent weak potentials is necessary, but not sufficient, to 
cause multi-angle decoherence. The negative sign in the 
Lorentz metric plays a crucial role because it determines 
that the flux term grows exponentially rather than oscil- 
lating. The occurrence of kinematical decoherence is a 
subtle feature of the nonlinear equations. 

We have studied the simplest possible case where deco- 
herence effects obtain. A more realistic case is provided 
by neutrinos streaming off a supernova core. In this case 
the fluxes of neutrinos and antineutrinos are different and 
geometry implies that the neutrino-neutrino interaction 
strength declines with the fourth power of radius. The 
evolution of an asymmetric system is far more compli- 
cated, even in the isotropic case, where the simple pen- 
dular motion turns into one of a spinning top that can 
precess, nutate, or simply swing, depending on the vv 
asymmetry. In the supernova context, transitions be- 
tween these modes occur as /i declines as a function of 
radius. Multi-angle decoherence cannot occur, for exam- 
ple, in regions where synchronized oscillations prevail. 

Numerical simulations suggest that for neutrinos 
streaming off a supernova core, the collective behavior 
characteristic of an isotropic gas prevails in many situa- 
tions, i.e., kinematical decoherence occurs, but is not the 
dominating feature of the system, contrary to what we 
have found here for the symmetric vv gas. Evidently it 
is not straightforward to apply the insights gained from 
our present study to the supernova case. 

We imagine, however, that the methods developed here 
may take us one step further to understand the supernova 
problem. We have formulated the equations of motion in 
a novel form adapted to the problem of multi-angle prop- 
agation, i.e., we have used a multipole expansion rather 
than the usual momentum modes and we have used the 
particle number and lepton number polarization vectors 
(or density matrices) as the fundamental variables. In 
this way, the crucial exponentially growing quantity, our 
Di, could be isolated as the primary cause of kinemat- 
ical decoherence. This approach may prove useful for 
understanding the supernova problem as well. 

For numerical three-flavor studies we have provided 
an expression for the conserved energy of the system in 
terms of the density matrices. The conserved energy and 
its three components E u , Eq and E\ provide a useful 
diagnostic tool for the behavior of the system. 

Even simple nonlinear systems can exhibit a surpris- 
ingly rich phenomenology. It is fascinating that the neu- 
trinos streaming off a supernova core provide an intrigu- 
ing case in point where many aspects of its behavior re- 
main to be understood. 
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